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^ , Abstract 



We present formulae for the momentum (z) distributions of D-mesons produced in neu- 
trino deep-inelastic scattering off strange partons. The expressions are derived within 
the variable flavor scheme of Aivazis et al. (ACOT-scheme), which is extended from its 
fully inclusive formulation to one-hadron inclusive leptoproduction. The dependence of 
the results on the assumed strange quark mass rris is investigated and the rris —>■ limit 
is compared to the corresponding MS results. The importance of 0{as) quark-initiated 
corrections is demonstrated for the rris = case. 



The momentuin {z) distributions of D-mesons from the fragmentation of charm quarks 
produced in neutrino deep-inelastic scattering (DIS) have been used recently to determine 
the strange quark distribution of the nucleon s{x, Q"^) at leading order (LO) [|T| and 
next-to-leading order (NLO) 0. A proper QCD calculation of this quantity requires 
the convolution of a perturbative hard scattering charm production cross section with a 
nonperturbative c ^ D fragmentation function Dc{z) leading at 0{as) to the breaking of 
factorization in Bjorken-x and z as is well known for light quarks [^. So far experimental 
analyses have assumed a factorized cross section even at NLO 0. This shortcoming has 
been pointed out in [Q and the hard scattering convolution kernels needed for a correct and 
complete NLO analysis have been calculated in the MS scheme with three massless flavors 
{u, d, s) using dimensional regularization. In the experimental NLO analysis in the 
variable flavor scheme (VFS) of Aivazis, Collins, Olness and Tung (ACOT) |Q for heavy 
flavor leptoproduction has been utilized. In this formalism one considers, in addition to 
the quark scattering (QS) process, e.g. W^s -^ c, the contribution from the gluon fusion 
(GF) process W^g -^ cs with its full r/is-dependence. The collinear logarithm which is 
already contained in the renormalized s(x, Q"^) is subtracted off numerically. The quark- 
initiated contributions from the subprocess W^s — > eg (together with virtual corrections) 
which were included in the complete NLO (MS) analysis in |^ are usually neglected in 
the ACOT-formalism. The ACOT-formalism has been formulated explicitly only for 
fully inclusive leptoproduction [0. It is the main purpose of this article to fill the gap 
and provide the expressions needed for a correct calculation of one-hadron (D-meson) 
inclusive leptoproduction also in this formalism. 

In the following we will stick closely to the ACOT formalism as formulated in |^ 
except that we are not working in the helicity basis but prefer the standard tensor basis 
implying the usual structure functions -Fi=i,2,3- We are not considering kinematical effects 
arising from an initial state quark mass in the W^s —> c quark scattering contribution, 
i.e., s(x,(5^) represents massless initial state strange quarks. This latter choice must be 
consistently kept in the subtraction term ^ to be identified below from the rris —* limit 
of the W~^g — > cs gluon fusion contribution. The fully massive partonic matrix elements 



have been calculated for the general boson-gluon-fusion process Bg -^ Q1Q2 in where 
B = 7*, W^, Z. When they are convoluted with a nonperturbative gluon distribution 
f7(x,yU^) and a fragmentation function Dq^{z), one obtains the GF part of the hadronic 
structure function Fi{x,z,Q'^) describing the momentum {z) distribution of a hadron H 
containing the heavy quark Q2'- 
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with the fractional momentum variables z = Ph-Pn/q-Pn and ( = pq^ -Pn/q-Pn, Pn and q 
being the momentum of the nucleon and the the virtual boson, respectively. The structure 
functions Fi{x,z,Q'^) generalize the usual fully inclusive structure functions Fi^XjQ"^), if 
one considers one-hadron (H) inclusive leptoproduction. The partonic structure functions 
fiix',C,Q^) are given by 
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where Am" = m2 — wi" , mi 2 being the mass of the heavy quark Qi^2- The kinematical 
boundaries of phase space in the convolutions in eq. (1) are 
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and axialvector {A) couphngs of the 7^(V^ — ^75) quark current enter via g± = V^ ± A^, 
Rq = V A . If the partonic structure functions in eq. (2) are integrated over C, the well 
known inclusive structure functions ||, |^ for heavy flavor production are recovered: 

SminK:^ ) 

where the fi{x', Q^) can be found in 
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In the following we will consider the special case of charged current charm production, 
i.e., mi = nisi m2 = f^tc {(1± = 2,0; i?g = 1 assuming a vanishing Cabibbo angle). Of 
course, all formulae below can be trivially adjusted to the general case of eqs. (1,2). The 
rris — > limit of the partonic structure functions in eq. (2) is obtained by keeping terms 
up to 0(1711) in the Ai, Ci and in C^max due to the singularity of the phase space integration 
stemming from (^ — > 1. One obtains 
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where PJIl^{x') = ^[a;'' + (1 - x')\ X = Q^/{Q^ + ml), Ci = 1/2, cs = x'/\ ,c, = 1 and 
the H? are the same as the dimensionally regularized MS {mg = 0) gluonic coefficient 
functions obtained in |^. The q arise from different normahzations of the fi and the Hf 
and are such that the infrared-safe subtracted [see below eq.(7)] convolutions in eq. (1) 
converge towards the corresponding ones in 0] as rus — *> if one realizes that x'/X = ^', 
x/\ = ^ = x{l + ni'l/Q'^). Taking also the limit mc — *> in eq. (5) gives -besides the 
collinear logarithm already present in eq. (5)- finite expressions which agree 0] with the 
massless results of 0. 

In the ACOT formalism the GF convolutions in eq. (1) coexist with the Born level 
quark scattering contributions F^ {x,z,Q'^) = ki s{^,^'^) Dc{z), ki=i^2,3 = 1, 2^, 2. The 
overlap between the QS and the GF contributions is removed by introducing a subtraction 
term (SUB) [H which is obtained from the massless limit in eq. (5) 



"^(a*^) 1 Ai 



.SUB _ j^^ ^^sKf^ In 4 
27r mi 



1 ^^' 



'"^Ki-'-i^'^pS''^ 



^, - DAz) . (6) 

The complete 0{as) structure functions for the z distribution of charmed hadrons (i.e., 
dominantly D-mesons) produced in charged current DIS are then given in the ACOT 
formalism [^ by 

pACOT ^ pQS _ pSUB _|_ pGF /^\ 

In Fig. 1 we show the structure function F^'"'^'^ at experimentally relevant ^ values 
of X and Q^ for several finite choices of vHs together with the asymptotic rdg — > limit. 
For Dc we use a Peterson fragmentation function |[Tl[ 

D,{z) = n[z{\- Z-' - ej{l - z)fY' (8) 

with e^ = 0.06 ^ 0] normalized to j^ dzD^{z) = 1 and we employ the GRV94(HO) 



parton distributions [|I4| with m,c = 1.5 GeV. Our choice of the factorization scale is 
fi^ = Q"^ + ml which ensures that there is no large \n.{Q'^ + ml})/ ^'^ present in the difference 
GF-SUB. As can be seen from Fig. 1 the effects of a finite strange mass are small and 
converge rapidly towards the massless MS limit provided m^ < 200 MeV as is usually 
assumed 0. 



In Fig. 2 we show the effects of adding the quark-initiated 0{as) correction from 
the process W~^s -^ eg (together with virtual corrections) to the asymptotic {rris —> 0) 
F^^^^, employing the CTEQ4(MS) densities fTD with m^ = 1.6 GeV. The 0{as) quark 
contribution is usually neglected in the ACOT formalism since it is assumed to be effec- 
tively suppressed by one order of a^ with respect to the gluon fusion contribution due 
to s{x, jj?)/ g{x, jj?) = 0{as). To check this assumption for the quantity F2{x,z,Q'^) we 
show, besides the full result, the contributions from the different subprocesses (using again 
/i^ = Q^ +m1). The W^g -^ cs contribution corresponds to GF-SUB in eq. (7). The 
quark-initiated 0{as) contribution has been calculated in the MS scheme according to 0] 
which is consistent with the asymptotic gluon-initiated correction in the ACOT scheme 
due to eq. (5). It can be seen that the quark-initiated correction is comparable in size to 
the gluon-initiated correction around the maximum of F2 . Since most of the experimen- 
tally measured D- mesons originate from this region the 0{as) quark contributions should 
not be neglected in a complete NLO calculation. It would be worthwhile to calculate 
these diagrams also in the ACOT scheme to study finite mg effects as has been done in 
this article for the 0{as) gluon contributions. 

To summarize we have given formulae which extend the ACOT scheme for the lep- 
toproduction of heavy quarks from its fully inclusive formulation to one-hadron inclusive 
leptoproduction. We have applied this formulation to D-meson production in charged 
current DIS and studied finite rris corrections to the asymptotic rris — > limit. The cor- 
rections turned out to be small for reasonable choices of rris ^ 200 MeV and we have 
shown that the m^ — ^ limit reproduces the dimensionally regularized MS (m^ = 0) 
gluonic coefficient functions |^. Furthermore we have investigated the quark-initiated 
0{as) corrections for rus = using the relevant MS fermionic coefficient functions |^. 
The latter corrections turned out to be numerically important at experimentally relevant 
values of x and Q^ ^ and should be included in a complete NLO calculation of charged 
current leptoproduction of D-mesons. 
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Figure Captions 

Fig. 1 The structure function F^'-^'-^'^{x, z, Q"^) as defined in eq. (7) using the GRV94(HO) 
parton densities |0 with irtc = 1.5 GeV and a Peterson fragmentation function [|ri|] 



with Ec = 0.06. Several finite choices for rris are shown as well as the asymptotic 
rUs — > limit. 

Fig. 2 The structure function F2{x,z,Q'^) for charged current leptoproduction of D- 

mesons at 0{as) for rris = using the CTEQ4(MS) parton distributions |T3[ and a 



Peterson fragmentation function [11| with Ec = 0.06. The full 0{as) result is shown 
as well as the individual contributions from the distinct quark- and gluon-initiated 
processes. 
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